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^ Abstract 

The equations for variation of the Stokes parameters and intensity of photons 
propagating in a medium, whose optical properties may be described by the permit- 
tivity tensor, are derived. Classification of different cases of photon propagation is 
Ph ■ suggested on a basis of these equations. 

in 

^ : 1 Introduction 

It is well known that polarized photons propagating in an optically active medium change 
initial polarization state 0,0). It has important meaning in a lot of fields of physics. We 
have pointed to some existing problems: 

a) propagation effects in a magnetosphere of pulsars || £| ; 

b) results of the circular polarization measurements of extragalactic radio sources || 
o and || (in our Galaxy); 

^ ■ c) cosmic birefringence (see and literature therein). 

d) problems of propagation x-rays |§ and high energy 7-quanta || [H| in single crystals 
and creation of the polarimeters for hard photon beams; 



O 



e) experimental observation of the birefringence in an electromagnetic vacuum |TT| , |12 

f) investigations of propagation of 7-quanta in laser waves |13|, 14] with the aim of 
■ polarization control on future 77-colliders. 

Besides, the possibility of propagation polarization phenomena of the radiation emitted 
from quark-gluon plasma must not be ruled out. 

Various phenomena (birefringence, Faraday rotation and etc.) of the visible light are 
well known and described in literature [p], 0- In @ the detailed consideration of these 
processes was carry out in anisotropic and gyrotropic media. However, the description of 
propagation was realized on a basis of the polarization states, whereas representation in 
terms of density matrices is more adequate in some cases. A parallel with well known ones, 
the new cases of propagation of photons were investigated in recent papers [[13], [L4], |L3| . In 



this connection the problem of classification of these cases and its completeness arises. 

On the other hand, the description of the optical properties of a medium with the help 
of the permittivity tensor is the universal and traditional method. This method may be 
used in wide energy ranges of propagating photons. The proof of validity of the method for 
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7-quanta and samples of its using one can find in [|16|, [L7|, pl| . Using of permittivity tensor 
presumes the linear response of media on perturbations. It is a very good assumption for 
many cases (however, see fL9| ). 



The aims of this paper are: firstly, the derivation of equations, describing all the cases of 
propagation of photons in media, whose optical properties may be described with the help 
of the permittivity tensor and, secondly, creation on a basis of these equations classification 
of the various cases of photon propagation. Our consideration is based on some suggestions 
such as smallness of the permittivity tensor components in compare to unit and straight line 
motion of the photons. It is not the most general case, of course, but the most interesting 
one for practical applications. 

It should be noted, that the problems of propagation of the photon beam were consid- 

H 



ered in some papers [20, 21 



the equations of the general type 



[22| . For example, in 

for propagation process were obtained. However, it seems that their using in specific ap- 
plications is difficulty, because the physical nature of the parameters in the equations is 
unknown. 



2 Equations 

We write the equations of the electromagnetic field in a medium in the following form Jl]. Q : 

rotB = -— — , div~D = , 

c dt 

1 <9B 

rotE = — , divB = 0, (1) 

c at 

where E is the intensity of electric field and D and B are the electric and magnetic induction 
vectors, t is the time, c is the speed of light. All the properties of the medium are reflected 
in the relation between B, E and D. Eqs.(|) would suffice to describe the propagation of 
photons in a medium and such a property as the intensity of magnetic field is not needed 
|2|. We represent the relation between D and E in the form 

Di{u) =e ij E j (u), {i,j = 1,2,3), (2) 

where = e'^ + ie'-j is the complex permittivity tensor and uj is the frequency of the 
photon beam propagating in the medium. 

We consider here only monochromatic photon beam, which propagates in defined di- 
rection. In this case conveniently to use in calculations the inverse t] a p of the tensor e^. In 
Cartesian coordinate system with one axis oriented along the direction of photon motion 
this tensor can be represented as two-dimensial tensor (a = 1, 2, /5 = 1,2) because of per- 
pendicularity of D to the photon wave vector. The obtained below equations are written 
in this coordinate system. 

In our consideration we assume that the photon beam move along straight line. It 
means that in transverse directions (relative to motion of the beam) a medium is uniform 
(0). In the most of practical cases it is true with very high accuracy. 

Knowing the tensor r] aj 3, one can find the refractive indices of propagating photons 
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n~ 2 



(Vu + %2) /2 ± \] (Vn - ??22) 2 / 4 + W72i , (3) 
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Thus, in the general case the photon beam propagates through the medium as the super- 
position of two electromagnetic waves with different refractive indices. In the general case 
the refractive indices are complex values. 

One can describe the polarization state either (of the two) wave by the use of the Stokes 
parameters X i; Y$ {i — 1 — 3) (AV values correspond to one wave and Y^ to another). These 



parameters are determined by the following relations [jM, IS 



= (4) 

X„Y, = (5) 

1 + KK* 

X 3 ,Y 3 = ^±. (6) 

1 + KK 

where k is the ratio of components of electric induction vector T)(Di, D 2 , 0) : 

D 2 mi \D 2 \ 

where k calculated with refractive index rii, (n 2 ) corresponds to Xj (Y^) -parameters (denote 
the K-values correspondingly as k% and k 2 ) and 5 is the phase shift between Di and D 2 . 

This ratio k can be reduced to zero or to the form k = ip (since |Di||Z?2| sin 5 = b\b 2 , 
where b\ and b 2 are the semiaxes of the ellipse and \p\ = b\jb 2 |]) by the rotation of the 
coordinate system around the photon wave vector (the third axis is constantly aligned with 
the wave vector). The first case corresponds to the propagation of a linearly polarized wave 
and the second case corresponds to an elliptically polarized wave; in addition, p > 0(p < 0) 
corresponds to left (right) - hand polarization of 7-quanta. 

From our consideration follows the important conclusion that these two waves (named 
as electromagnetic normal waves 0) with determined refractive indices and polarization 
states are the eigenfunctions of problem and any monochromatic electromagnetic wave 
propagating in the medium in defined direction is superposition of these two waves. 

One can obtain the following relations: 

KlK 2 = , (8) 

V21 

In the case of symmetric T^-tensor, the product k>ik 2 = —1 and X\ = —Yi, X 2 = Y 2 , X3 = 
— Y3. Another important case take a place when 7712 + 7721 = (but rju 7^ 0) and then 
kik 2 = 1, X\ = —Yi, X 2 = Y 2 , X 3 = —Y 3 . However, in general the absolute values of 
circular and linear polarizations for both the normal waves are not equal in between (i.e. 
\X 2 \^\Y 2 \ } Xl + Xl^Y? + Yi). 

Now we consider the case of photon propagation in the medium, whose optical prop- 
erties described by the tensor r\ a p at condition r\\ 2 + r\ 2 \ = 0, (7712 7^ 0). This case was 
investigated in [Oj] . The variations of intensity and Stokes parameters one can write in the 



following form: 

J y (x) = Ji(x) + J 2 {x) + 2J 3 (x) , (9) 



3 



d(x) = (X 1 J 1 (x) + Y 1 J 2 (x) +p 1 J 3 {x))/J 1 (x) , (10) 

= (X 2 J 1 (x) + Y 2 J 2 (x) +pMx))/J J (x) , (11) 

Ux) = (XMx) + Y 3 J 2 (x) +pMx))/J 1 {x) , (12) 

2 2X 3 2X 2 . . 

Pl = x~' P2 = ~^T' P3 = "xT' ^ ' 

where J 7 , ^1,^2^3 are the intensity and Stokes parameters of the photon beam on the 
distance equal to x. The partial intensities J{(x), (i = 1 — 4) have the following form: 

J x {x) = Ji(0)exp(-2Im(n 1 )cux/c) , (14) 

J2{x) = ^2(0) exp(— 2 lm(n 2 )ujx / c) , (15) 

J 3 (a;) = exp(— Im(ni + n 2 )ux/c){J 3 (0) cos(Re(ni — n 2 )ux/c) — (16) 

J4(0) sin(Re(ni — h 2 )ujx/c)} , 

Ja(x) = exp(— Im(ni + n 2 )ujx/c){ J%{ti) sin(Re(ni — h 2 )ujx/c) + (17) 

J 4 (0) cos(Re(rai — n 2 )ujx/c)} . 

The initial partial intensities are defined from the following relations: 

w n x 1-^16(0) +^26(0) +^36(0) (m , 

Jl(0) = 2(Xf + Xj) ' (18) 

x/ n x i-*ifr(Q)-*26(0)-x 3 &(0) / 1q n 

J2(0) " 2(Xf + Xf ) ' (19) 

J3(0) " 2(X| + Xf) ' (20) 

x 1 (x 2 e 3 (o)-X36(o)) , 91 * 

J4(0) 2(Xf + Xf) • (21) 

The relations between Xj and Yi values were used, because of this the Fj-values are 
absent in Eqs.([l8|)-(|2T|). Besides, we assume that J 7 (0) = 1. Now one can obtain the 



differential form of these equations: 

c dJ r 



--J (-S-Cfc-Bfc), (22) 



uj dx 

= -CC1C2 - - ^2 - 2?6, (23) 

a; ax 



C d£ 2 
C d£ 



C(l-ef)-B66 + ^i, (24) 
+ (25) 



3 _ an c2 

uj dx 



where 



A={r}' xl -r ] ' 22 )/2 } (26) 

B = (»&-»&)/2, (27) 

C = W2 - ^i)/2, (28) 

2> = (»& - (29) 

-(< + <) = fau + »&)/2. (30) 
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It should be noted, that these equations are true only at \5 a p — r\ a p\ <C 1 , where 5 a p is the 
Kroneker ^-function. This is due to the fact that the intensity of photon beam is calculated 
approximately (J ~ DD*). 

The differential form of equations for the general case of T^-tensor may be obtained 
by the similar direct calculations. However, another method is more conveniently to use 
for it. Knowing the relations for Stokes parameters in an anisotropic medium [|10| , [15], one 
can find the differential equations for this case. They have a similar form as Eqs.(22)-(25). 
Combining these equations and Eqs.(22)-(25) we get 

J 1 (-g-J 7 Zi-C&-B£ 3 ), (31) 

C&2 - B^ 3 - M2 - X>&, (32) 

^16-^26 + ^1-^3, (33) 

Hits-Ctes + Vh + Sb, (34) 
where 



C dJ ■y 



us dx 
c d£ 2 
uj dx 
c d£ 3 

uj dx 



uj 


dx 


= F{1 


~£ 


= C(1 


-6 2 


= B{1 





£ = W 3 + r4 a )/2, (35) 
F=(Vi2 + V2i)/2, (36) 



and A, B, C, V, Q are the same as in Eqs. (|26f) -(|30D. 

It is clear that these equations are required and they describe the propagation process 
in the general case. It obvious from following: 

1) the form of relations for intensity and Stokes parameters for the general case is similar 
to Eqs.(9)-(21). The main difference is contained that Ji(0), ^(O), ^(O), J 4 (0)-values have 
another mathematical representation. However, it is significant that these values are linear 
functions of Stokes parameters; 

2) the derivatives of Stokes parameters in Eqs.(32)-(34) are equal to zero when £1 = 
Xi, Y 1: £ 2 = X 2l Y 2 and £3 = X 3 , Y 3 . This fact supports the choice of A, T- parameters 
(see Eqs.(26)-(29) and (35)-(36)); 



3) Eq.(31) for intensity was obtained in |14| for the general case. 



Note that (— Q — J 7 ^ — — B£ 3 ) > for any real medium. 

The obtained here equations describe the variations of intensity and Stokes parameters 
for photon beam propagating in the medium, whose optical properties may be represented 
by the use of 7/ a/ 3-tensor. Besides, one can get the following equation for variations of 



polarization P = J£( + £2 + £3 of propagating photon 

c dP 2 



2uj dx 



l-P 2 )(-^i-C6-^ 3 )- (37) 
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3 Classification 



Obtained in previous section equations describe all the cases of propagation of photons in 
a medium. The difference between various cases is determined by the kind of r^-tensor. 
Here we consider briefly all the possible cases of the photon beam propagation in media. 

First and foremost it should be noted that the maximal number of values, which de- 
terminate the variations of photon state at its propagation in a medium, is equal to seven. 
However, the (/-value is responsible only for the intensity of photon beam (see Eqs.(|3"T])- 
([34]) ) . By this means only six (or less) values (.4, is required for description of the 

variations of photon polarization state [JS] • On the other hand, one can make the 8- value 
(or jF-value) equal to zero by rotation of the coordinate system around the wave vector 
of photon. It is simplify Eqs.(|3l])-(|34"|), but does not reduce the number of parameters (in 
this case the angle of rotation is such a parameter). For further consideration we select 
coordinate system in which £ = 0. 

It is obviously that medium is transparent if ^ = 0. This condition take a place 
at G = 0,.F = 0,C = 0,£> = 0. The excellent example of transparent medium is the 
monochromatic laser wave of not high intensity. The beam of 7-quanta propagates in 
its practically without interactions if the energy of the beam is less then the threshold 
of electron-positron pair production. Taking into account the importance of this case for 
practice, we write the solution of Eqs. (|3T]) -(|34D for transparent medium in the following 
form: 

Ux) = 6(0) cos0 + (6(0)X 3 - 6(0)X 2 ) sin0, (38) 
Ux) = -6(0)X 3 sin0 + X 2 £ 2 (0) + X 2 X 3 £ 3 (0) + (6(0)X 3 2 - £ 3 (0)X 2 X 3 ) C os0, (39) 
6(0) = 6(0)X 2 sin0 + X 2 X 3 £ 2 (0) + X 3 2 6(0) + (£ 3 (0)X| - £ 2 (0)X 2 X 3 ) cos0, (40) 

where = Re(ra 1 — n 2 )ujx/c and X 2 , X 3 are calculated with the help of Eqs.(f|)-(|6|). Specifi- 
cally, these equations describe the well-known cases of the light propagation as birefringence 
(X 2 = 0, X 3 = ±1) and Faraday rotation (X 2 = ±1, X 3 = 0). It easy to see that P = const 
in the transparent medium. Note that the absolute transparent medium is only physical 
abstraction 0. 

In general any tensor of second rank may be decomposed on the two parts: symmetric 
and antisymmetric. The antisymmetric part is dual to some axial vector, named as the 
vector of gyration. When the vector of gyration is equal to zero, the medium is named as 
anisotropic one. It means that the optical properties of anisotropic medium are determined 
by a symmetric //^-tensor, or, in other words, in the case under consideration the following 
relations are true: C = 0, T> = 0. We remind that S = due to the choice of the 
coordinate system. The anisotropic media may be classified by the help of the two types 
of ?7 Q/ g-tensor. In the first case T = 0, and JF ^ in the second case. The linearly 
polarized monochromatic laser wave is a simple example of the anisotropic medium of 
the first type ||14| . In this case the normal electromagnetic waves are linearly polarized 
along the principal axis of the ?7 Q ^-tensor and result of propagation of the initially linearly 
polarized photon beam in this medium is the effect of birefringence. This effect applies to 
transparent 3 = and absorbed B 7^ media. However, in absorbed media of considered 
type the initially unpolarized propagating photon beam becomes linearly polarized (at 
x — > 00). 
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q 


a n 




A ft 


a 3 — ±1, r 3 — — a 3 
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a n 

A-2, U 




A ft T 


V "ir^nv Y V Y 

1 2 — -^2 T U ? Y l — — A l? *3 — —A3 


(J 


G T 




V 


A^2 — -^-1? Y 2 — ^"2 


6 


G, D 




c,v 


A2 = ±1, Y 2 = —X 2 


7 


G + Ax,T 




A,V 


X 1 = Y 1 = 0,Y 2 = -X 2 ^ 0,F 3 = -X 3 


8 


G + A U D 




A,B,C,V 


Y 2 = -X 2 ^0,Y 1 = X 1 ,Y 3 = -X 3 


9 


G + A 2 ,D 




A,B,C,V,F 


\X 2 \?\Y 2 \,X* + XZ?Y? + Y? 



Table 1: The different cases of propagation of photons in media. 

The first letter(s) in the second column is the name of case, the letter after comma denotes 
the transparent (T) or absorbent (D) medium. Here I is an isotropic medium, A±, A 2 are 
anisotropic media of the first and second type, G is the gyrotropic medium. £ = due to 
the choice of the coordinate system. Symbol ± means that one of the two signs is used. 



The case of anisotropic medium of the second type take a place only for absorbed 
media. This case was considered in different media like single crystals ||10|| , plasma ||, 



laser wave |15] . In this case the normal electromagnetic waves are elliptically polarized. It 
is interesting that X 2 = Y 2 for this waves. The initially unpolarized propagating in this 
medium photon beam becomes linearly and circularly polarized. The simplest sample of 
the anisotropic medium of the second type for high energy photons is a dichromatic laser 
wave. (The dichromatic wave is a superposition of the two linearly polarized laser waves 
with different frequencies moving in the same direction and, broadly speaking, nonzero 
angle direction of polarization of these wave.) 

The case of isotropic medium is described by the following relations: A = 0, B = 
0, C = 0, T> = 0, T = 0. It is obviously that the initial polarization state of photon beam 
is conserved at its propagation in the isotropic medium. 

Now we consider the medium with the nonzero gyration vector. There are five cases 
of propagation of the photon beams. The first and second cases are the pure (transparent 
and absorbent, correspondingly) gyrotropic medium: ^4 = 0, £> = 0, 8 — 0, JF = and 
one or both of the values C, T> are not equal to zero. The normal electromagnetic waves 
are circularly polarized [X 2 = —Y 2 ). The Faraday rotation take a place if C — 0. The 
circularly polarized monochromatic wave is the sample of such a medium (C = 0, D / or 
both the values are not equal to zero) for propagating 7-beam with energy below (above) 
of threshold of pair production. 

The third and fourth cases are the gyrotropic medium with anisotropy of the first type 
{B = 0, C = 0, T = for transparent and T = for absorbent media). These cases are 
described with the help of Eqs.(|38|)-(|40|) and Eqs.(P)-([T2|) for transparent and absorbent 
media, correspondingly. The elliptically polarized monochromatic laser wave is the sample 
of such a medium. 

The last case is the gyrotropic medium with the anisotropy of the second type. The 
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normal electromagnetic waves for this medium have the different absolute values of circular 
and linear polarizations (I-X2I 7^ 1 1^1,^1 + X 2, 7^ Yj 2 + F 3 2 ). The sample of this medium 
is bichromatic combination of the two elliptically polarized waves moving in the same 
direction with the arbitrary angle between the semiaxis of polarization ellipses of these 
laser waves ( for propagating 7-beam above the threshold of pair production). 

The considered in this section cases of propagation are represented in table 1. In 
this table we assume conventionally that medium is transparent (T) independently of the 
quantity of ^-value. It allows to avoid the consideration of additional unimportant cases. 

Note, that the problems of C, P, T symmetries at propagation of the photon beam were 
investigated in EI . 



4 Discussion 

Eqs.(|3ll)-(|34]) describe all the cases of propagation of photons in uniform (the components 
of r^-tensor are constant along the direction of photon beam motion at fixed u) media. 
These differential equations of the first order have enough a simple form and dependent 
on differences of the corresponding components of r^-tensor. Besides, the equations were 
obtained at condition that \rj a p — 8 a p\ <C 1. One can find the some solutions of Eqs.(|3TD- 
(P^) depending on a kind of 7/ a( g-tensor |T^, |TJ]] . We do not know universal solution of these 
equations (the first integrals) describing all the cases of propagation. 

Eqs.(|3TD-(j34|) one can use for calculations of the propagation of photon beams in not 
uniform media, when the components of the ?7 Q/ g-tensor are functions of coordinate along 
direction of the photon wave vector. The condition, when this approximation is warrant, 
is smallness of dimensionless parameters c/(Ruj) relative to unit, where R is the distance 
of constancy of tensor components. In other words, these components should be fixed on 
the distances of photon wave length. 

Note, that in paper |13| the process of propagation of the high energy 7-quanta in the 



monochromatic laser wave of arbitrary polarization was considered. For this purpose the 



well known scattering amplitudes for elastic scattering light by light [16 were used. As a 



result the differential equations similar in form to Eqs.(E2])-(P^) were obtained. Comparison 



of Eqs.(p^)-(|25|) and similar ones in |13| allows to make conclusion that in this case the 
components of r^-tensor are linear combinations of the invariant helicity amplitudes for 
the forward light by light scattering. 

The propagating in an medium photon beam can be represented as a superposition of 
two normal electromagnetic waves with the different refractive indices. Because of this, 
in absorbent media one normal wave is absorbed to greater extent than other and after 
propagation of some distance only this wave would then be left behind. 

On the other hand, every normal wave propagates in such a way that its polarization 
state is conserved. Setting d^i/dx = 0, d&fdx = 0, d^/dx = in Eqs. ([3ll) -([34|) we 
can find the Stokes parameters of these waves (see Eqs.(||)-(|6|)). However, in the case of 
transparent medium the equations have infinite set of solutions and the third equation is 

+ el + el = 1. 



Eq.(|37[) is the consequence of Eqs.(pT])-(p^). The total polarization \P\ is decreased or 



increased in the relation of sign of the right side of Eq. fl3T|) . In the case, when £>, C, T are 
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equal to zero P = const. On the other hand, it is easy to get for the absorbent medium: 



,2, s , a-n 2 )exp(25wi/c) 



pi{x) = , _ , U^==^ , (4!) 

where Pq = P(0) is the initial total polarization. Note, that J 7 (0) = 1. Besides, the value 
—Quj/c is the inverse mean path of the unpolaraized photons in the medium and may be 
measured in principle. The obtained simple relation is universal and may be useful for 
experimental determination of photon beam polarization. 

We have considered the propagation of monochromatic photons in media. The absorp- 
tion of photon means one of the two possibilities: a) scattering; b) vanishing. Compton 
effect and pair production are the samples of these possibilities for 7-quanta. As a result, 
the cascade process take a place in media. As a rule for description of cascade processes 
the kinetic equations are used p5 |. The kernels of these equations are the differential prob- 



abilities (in particular, per unit length) of the possible elementary processes. From this 
point of view Eqs.(]3]J)-(]34]). may be useful in similar calculations. 

We understand that our classification (see Table 1) is a matter of convention to a degree. 
For instance, cases 5 and 6 from table 1 at very small C-value give practically the same 
description of linearly polarized photon-beam propagation on short distances. However, 
even so, photon beam obtains a small value of circular polarization (in case 6). Note, 
that on infinity this beam becomes completely circularly polarized, whereas in case 5 beam 
remains always linearly polarized. 

Besides, case 4 at A = ±JF and B = may be considered as special one. Under this 
condition the refractive indices of normal waves are equal in between and X2 = Y2 (see for 
detail @). However, we interpret its as limit of case 4. 

It should be noted that the analytical solutions of Eqs.(PT|)-(P4]) are exist and their full 
number is less than the number of cases in Table 1. The universal solution for cases 1-3, 
5-7 one can found in the appendix of paper ||14||. The propagation in anisotropic medium 



of the second type (case 4) is considered in [pi], [15| . Eqs.(|9"D-(|T2D describe the propagation 



process in the gyrotropic medium with the anisotropy of the first type. 



5 Conclusions 

The derived equations for intensity and Stokes parameters of the monochromatic photons 
propagating in an uniform medium, whose optical properties may be described by the per- 
mittivity tensor, is a central result. The equations have the simple form and the variation 
of the Stokes parameters is determined only by the differences of corresponding compo- 
nents of the inverse r] a p of the permittivity tensor. It is shown that the equations may 
be used under certain conditions for calculations of the photon beam propagation in not 
uniform media. 

On the basic of these equations the classification of different cases of photon propagation 
are suggested. 

The problem of generalization of Eqs.fl3T|)-(|3"4|) for arbitrary quantities of permittivity 
tensor components is of interest and requires further consideration. 

Our results may be useful for investigations in astrophysics, and in some areas of high 
energy and nuclear physics, and physics of condensed media. 
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